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The dose rate in a rotationally- symmetric radiation field near a space power reactor 
is assumed to be inversely proportional to separation distance squared and directly propor- 
tional either to a constant or a constant times the sine-cubed of a polar angle. Constant 
velocity motion ( constant speed and direction) is examined in both cases in both two- and 
three-dimensional geometries. The two-dimensional geometry occurs when the line of 
motion and the field axis of symmetry are co-planar. The dose integral in the sine-cubed 
region may be integrated directly, but a more useful form is obtained after a change of vari- 
able. A coordinate system rotation greatly simplifies the results. The three-dimensional 
problem is integrated after a change of variable. Finally, tables of normalized functions are 
presented and discussed. 


INTRODUCTION 


£ = constant , 


A key problem associated with the study of 
nuclear energy sources in space is the evaluation of the 
leakage radiation reaching the astronauts or other radio- <j> 

sensitive detectors. This paper discusses a solution to 
the problem of constant velocity motion near a space and 

power reactor with a specific leakage radiation pattern. 

The term ’’constant velocity” is intended to be under- dp 

stood in the vector sense implying constant speed and d t 

constant direction. 


constant , 
constant , 


V - constant . 


The radiation field description was derived from 
material presented by Atomics International in Refer- 
ences 1 and 2 and is illustrated in Figures 1 and 2. 
Figure 3 depicts an isodose contour in the (x, +y) half- 
plane. The radiation source is represented as a point 
source located at the origin O. The radiation field is 
assumed to be comprised of three regions. Region I 
represents the field due to an unshielded isotropic 
source or leaking through a relatively thin shield. 

Region III represents the field behind a thick ’’biological 
shield” of half- angle e. In both regions, the shield is 
assumed to be of approximately uniform attenuation with 
azimuth angle so the dose rate is a function of distance 
alone. Region II represents the transition between 
regions I and HI. Because of variations in attenuation 
and/ or source strength, the dose rate has an azimuthal 
dependence. The field studied here varies as the sine- 
cubed of a polar angle (p measured from an abscissa 
inclined an angle £ to the x-axis. 


In region II, the dose rate is 

D - K 1 si " 3 £ 

n x* - 


but 


thus. 


0 - (a - £) ; 


a a sima _ sin a sin a* 

D u = A — p— + b -p- + c p 


+ d £2^: 


where 


A = Kj cos 36 , 


MOTION WITHIN A SINGLE PLANE 


B = Kj cos 6 sin 2 6 , 


Consider first the case where the detector moves 
with constant velocity along a line p contained within a 
nonrotating plane which also contains the x-axis. Let us 
first examine the case where all motion is in region II, 
or 


u - 6) < 0 s J . 

Let the motion be outward along p from an initial posi- 
tion Pj. Then, from Figure 4, 


C = Kj sin 6 cos 2 6 , 

D - Kj sin 36 , 



and 


r 


= R 


II 


r( P^ = i’j = constant 


Consider the first term. Referring back to 
Figure 4 and introducing initial conditions, 
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r sin a = r t sin a x + p sin £ . 

Also, 

r 2 = p 2 + rj 2 + 2 prj cos ( £ - 04 ) 


By integrating and inserting the limits, we have the com- 
plete first term, 

_ A 3(cos a i - cos a 2 ) -<(cos 3 a i - cos 3 a 9 ) 

IIA 3r t V _ sin (£ - a x ) 



{r A sin q 1 + p sin £) 3 dp 

[rj 2 + 2prj cos (£ - Oj) + p 2 ] ^ dt 

The limits for the complete integral along the line are: 

P(l) = 0 

and 

p( 2 ) = 00 . 

By integrating and inserting the limits and simplifying, 
the equation above eventually yields 

D = A I (sin 3 g t + sin 3 £) 

IIA 3r t V \ [1 + cos (£ - Qfj)] 

(sin a f + sin £) 3 I 
+ [1 + cos (£ - a t ) 2 ] I 


The second term of the expanded dose rate equa- 
tion is 

• _ B sin a 

IIB " r 2 

which yields immediately 

D _ B (cos a i - cos a 9 ) 

IIB r t V sin (£ - c^) 

using the limits above. 

The doses and may be evaluated using 

the arguments above but with the complementary angle 
a\ The complete dose expression can now be written as 

D II = 3rjV sin (f - Tj {cOS 36(3 (COS “l " c°s a,) 

- (cos 3 a i - cos 3 a 2 )] 


We now write 


IIA V 


CHANGE OF VARIABLE 


+ 3 cos 6 sin 2 6 (cos a cos a 2 ) 


This general form is somewhat cumbersome to - 3 sin 6 cos 2 6 (cos af - cosa 2 ') 

handle. A more useful, if somewhat restricted, form 

may be obtained through a change of variable before + sin 36 [3 (cos - cos aj) 

integration. Referring again to Figure 4, we observe 

that - (cos 3 a/ - cos 3 c^’)]} 


r cos a = r t cos a t + p cos £ . 


Taking this in conjunction with the expressions above, 
we find that 


r 


r sin ( £ - g t ) 
1 sin (£ - a) 


This iorm is valid except in the case of radial motion 
discussed below. An additional minor constraint is that 
> e . 


This expression can be greatly simplified by per- 
forming a coordinate system rotation through an angle a. 
Thus, 


sin (q - «<) 

P 1 sin (£ - a) 9 

and 

dp = sin ( £ - gj) da 
dt sin 2 (£ - a) dt * 

Substituting into the dose rate expression, 

• = sin 3 a da 

IIA sin (£ - ctj) dt 

The limits will be in general 

a(l) = 


and 


a * - a - 5 = <p , 

f* = £ - 6 , 

6 * = 0 , 


and, constraining the point Pj , to lie along the ax 
line, 


D n (0 ,<*>) 


K t cos 3 0-3 cos 0 + 2 
3rjV sin (£ - 6 ) 


= 6 


where the parentheses on the left are used to indicate the 
angular limits. The radial case may be integrated 
directly from the dose rate equation. It is unaffected by 
the coordinate system rotation; 


D U (0,<O = f 


J. 

r l 


1 

r 


sin 3 $ 


a{2) = a 2 . 




The case of constant velocity co-planar motion in 
an inverse square field has been solved before and is 
included here only for completeness and notational con- 
sistency. Referring to Figure 5, we again let the motion 
be along a line p , and, as before, 

r(P t ) = rj = constant , 

£ = constant , 

aj = constant , 

<t>\ = constant , 

and 

^ = constant = V . 
d t 

The dose rate is given by 




We observe that 

rj = f - K - a) . 

r = r* cos r)x sec rj , 
p ~ r t cos T] t (tan 17 - tan tjj) 


and 


FIGURE 5. CONSTANT VELOCITY COPLANAR MOTION 
IN AN INVERSE SQUARE FIELD 


ft = r i cos Tj, sec 2 tj fj ; 

thus, 

d = jl.t - . 

r t V [sin (£ - a t ) J 

MOTION IN THREE DIMENSIONS 

If isodose contours such as in Figure 3 are 
rotated through 2tt radians around the x-axis, a complex 
three-dimensional radiation field will be produced. 
Assume constant velocity motion along an arbitrarily 
oriented straight line, as in Figure 6 . Again, let us 
first examine the case of motion within region II along 
the line P A P = p outward from an initial position Pj. 
The radiation field polar angle is again 



<p = a - 6 , 

and, as before, 


FIGURE 6. CONSTANT VELOCITY MOTION 
IN THREE DIMENSIONS 


D n 


n 


= A 


sin 3 a 
r2 


+ B 


sin a 
r 2 


+ C 


sin a 1 sin 3 a' 
r 2 + r 2 


The direct integral of the cubic term has not been 
obtained; however, the term has been successfully inte- 
grated after performing a variable change similar to that 
in the restricted solution of the co-planar case. 
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The radius vector to the initial point can be 

written 


- .a ^ a 

r t = (l cos Oj + j cos + k cos y^rj , 

where r t = Jr^ and cos Uj, cos /3 1# and cos y t are 
direction cosines. Similarly, the vector traced by the 
moving point is 


p = (1 cos a + 
P 


j A cos fi 


£ cos y p )p 


Since the direction is constant, p/p is constant. The 
angle between these two vectors is also a constant, 


cos x = cos on cos a + cos 61 cos 8 
P Hl P p 



+ cos yj cos y 


Figure 7 ^illustrates the nonrotating plane contain- 
ing the vectors r, r lf and p. This plane contains the 
variable angle 0 between the vector r and the (x,y) 
plane. Consider the nonrotating plane containing the 
x-axis and and the rotating plane containing the 
x-axis and r. In the former, there is the constant angle 
and in the latter the variable angle a . If a unit 
sphere is constructed at the origin, these three angles 
become the sides of a spherical triangle. If the face 
angle opposite side a is called X, then 

cos a = cos cos (0 - 0j) 

+ sin a t sin (0 - 0 t ) cos X . 

But, X is also the angle between the planes containing 
POPj and PjOXj. Thus, 

cos Or - X) = i 

\ r iP / \ r i / sin X sin a, 

and 

cos a sin 2 - cosa t (cosp cosp t + cosy cosy^ 

cos X P P 

sin X sin a t 

Continuing, we simplify 

cos a = k cos [0 - (0j - j p)] , 

where the phase angle ip is defined by 
tan ip = tan cos X 
and the modulus k is found from 

k 2 = cos 2 a t + sin 2 a t cos 2 X . 


FIGURE 7. PLANE OF MOTION IN THREE-DIMENSIONAL CASE 


p = r si ? 

p 1 sin (£ - e) ’ 

and 

d£ _ sin (£ - dj ) d 0 
dt 1 sin 2 (£ - 0) dt 

The dose rate is now 

D = -A- U - k 2 cos 2 [0 - (0j + ip)]}^ d0 
UA v{V sin (£ - 0) dt * 

Let 

sin 2 fj, = cos 2 [0 - (0 t + ip)] 

and 

A = (1 - k 2 sin 2 n) ^ , k 2 < i 

and the expression integrates to 

d iia = sTTv sin a - e7) [k2A sin " cos ** 

- (1 - k 2 )F(k,n) 


We note in passing that $ lf ip, and k are defined by (2) 

initial conditions only and that k 2 < i. + 2 < 2 - k 2 )E(k,/x)]^ » 

Returning to Figure 7 , we may write where F(k,^) is the incomplete elliptic integral of the 

first kind and E(k,pi) is the incomplete elliptic integral 
r = r A sin il Z-gj) ^ of the second kind. The limits of integration are indi- 

sin (£ - 0) ’ cated outside the bracket on the right. 
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The complete expression for the three- 
dimensional case could now be constructed using the 
components above. Regrettably, no simplifying step 
such as the coordinate system rotation used in the planar 
case has been found. 

The restriction against radial motion applies 
here as in the planar case. However, as radial motion 
is necessarily co-planar with the x-axis, the special 
form given earlier will apply. 

APPLICATIONS 

Many problems involved with close-in operations 
near a nuclear power source can be adequately approxi- 
mated by assuming constant velocity co-planar motion. 
The integral dose in this case can most easily be found 
from the expression for D^CO, <t>) developed above. The 

expression can be normalized as 

( , rjV sin (£ - 6) (cos 3 <f> - 3 cos <p + _2) 

D U ( 9 (P> Kj 3 

This expression is tabulated as a function of the angle <j> 
in Table 1. This table is organized in standard math 
table format and allows five significant figure values of 
the function to be read corresponding to three significant 
figure values for the angle <f>. Linear interpolation is 
used for intermediate values. 


TABLE 1. (Concluded) 


0123456789 

45.0 07741 07803 07865 07928 07991 08054 08117 C8l8l 08245 08310 

46.0 08374 08440 08505 08571 08637 0«703 08770 08837 08905 08972 

47.0 09041 09109 09178 09247 09316 09386 09456 Q9527 09597 09669 

48.0 09740 09812 09884 09956 10029 10102 10176 10250 10324 10398 

49.0 10473 10549 10624 10700 10776 10853 1093Q 11007 11084 11162 

50.0 11241 11319 11398 1147B 11557 11637 11718 11798 11879 11961 

51.0 12043 12125 12207 12290 12373 12457 l254o 12625 12709 12794 

52*0 12879 12965 1305i 13137 13224 i33ll 13398 13486 13574 13662 

53.0 13751 13840 13929 14ol9 14109 14200 1429 q 14382 14473 14565 

54.0 14657 14750 14843 14936 15030 15124 15218 15313 15408 15503 

55.0 1559? 15695 1579g 15888 15986 16083 16181 16279 16378 16477 

56.0 16576 16676 16776 16876 16977 17078 l7l?9 17281 17383 17485 

57.0 17588 17691 17795 17898 18003 18107 18212 18317 18423 18529 

58.0 18635 18742 18849 18956 19o64 19172 i928q 19389 19498 19607 

59.0 19717 19827 19937 ?004a 20159 20271 20383 20495 20607 20720 

60*0 20833 20947 21061 2H75 21290 214o4 21520 21635 21751 21867 

61.0 21984 22101 22218 22336 22454 22572 22691 22810 22929 23049 

62*0 23169 23g89 234l 0 23531 23652 23773 23895 24ol8 24140 24263 

63.0 24387 24510 24634 24758 24883 25008 25133 25259 25385 255ll 

64.0 25638 25764 2589? 26019 26147 26275 26404 26533 26662 26791 

65.0 26921 27o5l 27181 $7312 27443 27575 27706 27*38 27970 28103 

66.0 28236 28369 28503 $8637 28771 28905 29040 29175 29310 29446 

67.0 29582 29718 29855 29992 30129 30266 30404 30542 30681 30819 

68.0 30958 31098 31237 31377 31517 31658 31798 31939 32081 32222 

69.0 32364 32506 32649 3279i 32934 33078 33221 33365 33509 33654 

70.0 33798 33943 34088 34234 34380 34526 34672 34*19 34966 35113 

71.0 35260 35408 35556 35704 35852 36001 36150 36299 36449 36599 

72.0 36749 36899 37049 37200 37351 37502 37654 37*06 37958 38110 

73.0 38263 38415 38568 38722 3887$ 39029 39183 39337 39491 39646 

74.0 39801 39956 4oU2 40267 40423 40579 40735 40*92 41049 41205 

75.0 41363 41520 4167* 41836 41994 43152 42310 42469 42628 42787 

76.0 42946 43106 43266 43*26 43586 43746 43907 44068 44228 44390 

77.0 44551 44713 44874 45Q36 45198 45361 45523 45686 45849 46012 

78.0 46175 46339 46502 46666 46830 46994 47158 47323 47487 47652 

79*0 47817 47983 48148 48313 48479 48645 488U 4«977 49143 493l0 

80.0 49476 49643 49Blo 49977 50144 50312 50479 50647 50815 50983 

81. 0 51151 51319 51487 51656 51825 51993 52162 52331 52500 52670 

82.0 52839 53o09 53178 53348 53518 53688 53858 54029 54199 54369 

83.0 54540 54711 54882 55f>53 55224 55395 55566 55737 55909 56080 

84.0 56252 56424 56595 56767 56939 571U 57284 57456 57628 57801 

85.0 57973 58i46 58318 58491 58664 58837 59olO 59i83 59356 59529 

86.0 59702 59876 60049 60222 60396 60569 60743 60917 61090 61264 

87.0 61438 61612 &l78fe 61960 62133 623 q7 62482 62656 62830 63004 

88.0 63178 63352 63527 63701 63875 64050 64224 64398 64573 64747 

89.0 64922 65o96 65271 65445 65620 66794 65969 66l43 66318 6649? 

0123456789 
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Values from this table give the integral dose 
along a line originating on the <p = 0 or - 6 line 
and proceeding outward until a - 6 = (p. The inclina- 
tion of the line is accounted for in the normalizing con- 
stants. Integral dose accrued between the points at 
and <j> 2 on the same line as [D^O, (p z ) - D^O - 4>i)]. 

This allows one to correct for the small differ- 
ence between a = e , the closest the axis may be 
approached in region II, and a = 6 , the abscissa of 
the extended field contours. 

Another useful form results from constraining P 4 
to lie along the <p = ^ line; thus, 

A 



This dose is related to the tabulated dose by 

D u(?’*) = [ D n( 0, f) " D n (0,< * >) ] ' 6) • 

The three-dimensional problem in region II con- 
tains too many variables to allow simple parametric 
representation. To obtain sample results, a specific 
field pattern was chosen and values of the integral above 
were obtained. The field axis offset angle, 6 = 8. 663 
degrees, corresponds to the field illustrated in 
Figures 1 and 2. Generalized values of the resulting 
integral are shown in Table 2 as functions of the modu- 
lus k and the field angle p defined above. The table 
is read in a manner similar to Table 1; thus, 

D n ( k .Mi.M2) = D n (k,/i 2 ) - D n (k,^i 1 ) , 

where it is understood that the normalizing factor 
[r t V sin (£ - $)/k] has been divided out of both sides 
of the expression. 


TABLE ii. NORMALIZED DOGE INTEGRALS FOR A THREE-DIMENSIONAL CASE, 
4 - 8. 667 DEGREES 
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